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摘 要 : 以 塔里木 胡杨 林 国 家 级 自然 保护 区 为 研究 对 象 , 基 于 植被 WPP 气象 .土地 利用 /覆盖 河流 
等 数据 ,采用 趋势 分 析 、Mann-Kendall 突变 检验 、 地 理 探测 器 等 方法 ,探究 了 塔里木 胡杨 林 国 家 级 自 
然 保 护 区 2000 一 2015 年 植被 NPP 时 空 变化 特征 及 其 影响 因素 。 结 果 表 明 :(1) 在 空间 尺度 上 , 塔 


5.16 ~303.87 gC © m”? 


里 木 胡 杨 林 国 家 级 自然 保护 区 2000 一 2015 年 年 均 植被 NPP 为 32.25 gC m°? 
“a 之 间 ; 年 均 NPP 呈现 出 以 塔里木 河 于 流 为 带 向 周边 波动 说 减 的 空间 


Ea] 在 


“a ,变化 范 


分 布 特征 。(2) 在 时 间 尺 度 上 ,16 a 间 保 护 区 植被 NPP 呈现 波动 增长 趋势 ,年 均 增加 值 为 0.523 8 


eC +m 


“a ,在 2001—2002 年 和 2011—2012 年 出 现 突变 性 上 升 ,2007 一 2008 年 出 现 突变 性 下 
降 。(3) 影响 植被 NPP 分 异 的 核心 因素 为 土地 利用 /和 窗 盖 、 燕 散发 、 降 水、 河流 缓冲 区 等 , 且 由 多 医 


子 协同 作用 造成 ;同时 ,地 下 水 埋 深 对 植被 NPP 变化 有 着 重要 影响 。 


关 键 词 : 


植被 净 初 级 生产 力 ( Net Primary Productivity ) 
是 指 在 自然 环境 下 ,绿色 植物 通过 光合 作用 在 单 
位 面积 、 单 位 时 间 内 积累 的 有 机 物 数 量 , 是 光合 作 
用 所 产生 的 有 机 质 总 量 扣除 自 养 呼吸 后 的 剩余 部 
A. NPP 作为 陆地 生态 系统 物质 循环 和 能 量 流 
动 的 基础 ,表征 陆地 生态 系统 生产 能 力 和 质量 状 
况 , 是 评价 陆地 生态 系统 健康 和 可 持续 发 展 、 估 算 
地 球 支 持 能 力 的 重要 指标 ”“ 。 随 着 全 球 变化 研 
究 的 不 断 深入 ,植被 净 初 级 生产 力 在 人 研究 全 球 变 
化 对 生态 系统 的 影响 、 响 应 和 对 策 中 ,成 为 一 项 不 
可 缺少 的 重要 指标 。 因 此 ,研究 植被 NPP 时 空 变 
化 状况 是 评价 生态 系统 健康 程度 及 稳定 性 的 重要 
内 容 ,为 全 球 性 或 区 域 性 的 生态 环境 问题 研究 提 
供 了 重要 依据 ” 。 

NPP 作为 反应 陆地 生态 系统 对 全 球 变 化 响应 
的 重要 指标 "| ,是 国内 外 学 者 长 期 关注 的 领域 之 
一 。 目 前 ,国内 外 学 者 的 研究 内 容 主 要 包括 两 方面 : 
一 是 植被 NPP 动态 演变 特征 ,人 研究 尺度 从 全 球 、 全 
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国 到 区 域 . 省 市 均 有 涉及 7 ;国内 学 者 的 研究 
区 域 主要 集中 在 长 江 黄河 流域 西南 地 区 南方 地 
区 ,东北 地 区 .内 蒙古 地 区 等 湿润 半 湿 润 地 区 … ， 
近年 来 干旱 半 干 旱 区 植被 NPP 动态 变化 也 受到 学 
者 '"“ 沾 的 关注 和 重视 ;二 是 植被 NPP 空间 分 异 的 
影响 机 理 研究 ,学 者 们 主要 选取 气温 .降水 等 影响 因 
素 , 运 用 相关 性 分 析 、 回 归 分 析 、SVD 分 析 等 方 
法 "9 ,分析 各 影响 因素 与 植被 NPP 的 相关 性 ， 
揭示 植被 NPP 空间 分 异 的 原因 。 纵 观 已 有 研究 发 
现 , 多 数 文献 基于 大 中 尺度 研究 植被 NPP 时 空 演变 
趋势 ,往往 忽略 了 保护 区 等 小 尺度 植被 NPP 变化 的 
细节 特征 与 规律 ;此 外 ,学 者 们 对 植被 NPP 的 影响 
机 理 虽 进 行 了 广泛 研究 ,但 主要 侧重 单个 或 少数 几 
个 影响 因素 ,缺少 对 影响 因素 重要 程度 及 交互 作用 
的 研究 。 

塔里木 胡 杨 林 国 家 级 自然 保护 区 是 世界 上 原始 
胡杨 林 生 长 最 好 、 保 存 最 完整 的 区 域 ,是 塔里木 河流 
域 重要 的 生态 屏障 ,对 于 维护 荒漠 脆弱 的 生态 平衡 
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具有 重要 的 作用 "| ,研究 全 球 气候 变化 背景 下 塔 里 
木 胡杨 林 国 家 级 自然 保护 区 植被 NPP 时 空 变化 及 
其 制约 因素 显得 尤为 重要 。 对 此 ,我 国学 者 开展 了 

一 系列 的 相关 研究 ,但 研究 大 多 从 保护 区 胡杨 径 向 

EKE 植被 覆盖 "胡杨 种 群 结构 与 格局 
等 方面 进行 分 析 , 对 于 保护 区 植被 NPP 空间 分 布 特 
征 及 其 影响 因素 并 没有 进行 详细 分 析 。 基 于 此 ,本 
研究 以 2000 一 2015 年 逐 月 植被 NPP 数据 为 基础 ， 
综合 应 用 趋势 分 析 、Mann-Kendall 突变 检验 以 及 空 
间 分 析 方 法 ， te A re 自然 保护 区 植 
被 NPP 时 空 演 变 特征 及 其 影响 因素 进行 系统 分 析 ， 
ee meres 参 断 及 保护 管理 成 效 评 
估 提 供 科 学 依据 。 


1 研究 区 概况 


塔里木 胡杨 林 国家 级 自然 保护 区 位 于 新 疆 巴 
Tr tS rh AA ALR 轮 台 县 、 库 尔 勒 市 境内 
(图 1) ,地 处 塔克拉玛干 沙漠 北 缘 ,位 于 我 国 最 长 
E FOR Tay AY Pie, BE E FLOR Hh AB 
段 。 保 护 区 东西 长 109.7 km ,南北 宽 47.1 km, 总 
面积 39. 54 x 10° hm’ ,是 干旱 荒漠 区 以 胡杨 林 为 主 
的 森林 生态 系统 类 型 自然 保护 区 ,属于 塔 克 拉 玛 
干戈 壁 荒漠 .大陆 性 荒漠 . 半 荒 漠 生 物 群 落 的 典型 
代表 。 保 护 区 属于 典型 的 温带 大 陆 性 平原 区 荒漠 
气候 ,区 内 温差 变化 大 ,多 年 平均 降水 量 为 45.2 
mm ,蒸发 量 却 高 达 1 887 ~2 910 mm, = 10 % 年 积 
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温 4 040 ~4 210 %C ,年 日 照 时 数 为 2 573.5 h ,无 
霜 期 218 do 


2 数据 来 源 与 研究 方法 


2.1 数据 来 源 

植被 NPP、 土 地 利用 / 禾 盖 数据 来 源 于 全 球 干 
旱 生 态 系 统 国 际 大 科学 计划 ( Global-DEP ) 提供 的 数 
据 产 品 ,该 产品 从 方法 上 把 遥感 数据 .生态 系统 监测 
数据 .野外 调查 数据 及 案例 研究 结合 起 来 , 逐 月 估算 
每 月 WPP ,空间 分 辩 率 为 230 m, 时 间 跨 度 为 2000 一 
2015 年 。 气 温和 降水 数据 来 源 于 中 国 科学 院 资源 
环境 科学 数据 中 心 (http: // www. resdc. en) ,该 数据 
采用 ANUSPIN 插值 ,插值 结果 比 协同 克 里 金 插 值 提 
高 了 13.23% , ERPS HE T 20.13% ; 
气温 和 降水 时 间 跨 度 为 2000 一 2015 年 ,时 间 分 辩 率 
为 年 ,空间 分 辩 率 为 1 km。 东 散发 数据 来 源 于 国家 
生态 系统 观测 研究 网 络 科技 资源 服务 系统 (http: / 
rs. cern. ac. cn/ ) ,该 数据 集 包 括 1 ~ 12 ARE RE 
散发 数据 ,空间 分 辨 率 为 1 km。 土 壤 类 型 数据 来 源 
于 寒 区 旱 区 科学 数据 中 心 (http : // westde. westgis. 
ac. cn/) ,根据 FAO-90 土壤 分 类 系统 ,对 土壤 类 型 
进行 重 分 类 。 植 被 类 型 参考 张 新 时 主编 的 中 华人 民 
共和 国 植被 图 (1:1 000 000) ,河流 等 基础 数据 来 源 
于 1:1 000 000 全 国 基础 地 理 数据 集 。 所 有 数据 统 
一 采用 Krasovsky_1940_Albers 投影 坐标 ,通过 重 采 
样 获取 空间 分 辨 率 为 250 m 栅 格 数据 。 
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图 1 研究 区 示意 图 


Fig.1 Location of the study area 


tió WN 


2.2 研究 方法 

2.2.1 趋势 分 析 法 ”以 象 元 为 基本 单位 ,采用 一 元 
线性 回归 方程 拟 合 2000—2015 年 共 16 a 的 塔里木 
胡杨 林 国 家 级 自然 保护 区 年 均 NPP 的 变化 趋势 , 反 
映 不 同年 份 植被 NPP 的 变化 趋势 及 特征 ,其 计算 公 
式 为 : 


式 中 :kj 为 单个 象 元 多 年 回归 的 斜率 ;n 为 年 数 ; 
NPP, 表示 某 一 点 象 元 第 i 年 的 NPP 值 ,其 中 i=1， 
2,3,…,16。 当 Kop >0 时 ,表示 随时 间 i ASM IN, 
NPP $ EFR; M ,<0 时 ,表示 随时 间 i Hy 
加 ,NPP 呈 下 降 趋势 ;1 1 的 值 越 大 ,表示 NPP 上 
升 或 下 降 的 速率 越 大 。 利 用 NPP 变化 率 表征 16 a 
间 VPP 的 变化 程度 ,其 计算 公式 为 : 


k 
NPP, DL 


me = gpp Xn X100% (2) 


INEP : NPP pange A AI RRKZR ; NPP mean BEAN 2000 一 
2015 年 共 16 a 年 平均 NPP;n 为 年 数 。 

2.2. 2 Mann-Kendall 趋势 分 析 及 突变 检验 
Mann-Kendall 突变 检验 常用 于 分 析 气 象 数 据 时 间 序 
列 趋势 ,采用 Mann-Kendall 突变 检验 分 析 塔 里 木 胡 
杨 林 国 家 级 自然 保护 区 近 16 a 年 均 NPP 的 变化 趋 
A ,并 分 析 其 突变 点 。 计 算 公 式 为 : 


k 
d, = 2 %; 


式 中 :n 为 样本 容量 ;x; 为 时 间 序 列 。 的 均值 和 方 
差分 别 为 : 


(2<k<n) (3) 


E(d,) =k(k -0. 25) (4) 
Var( d,) -FED (2E +5) (5) 


标准 化 后 得 到 UF: 


d,- E(d,) 
UF n) = 一 天 一 一 一 

VVarl d,) 

由 公式 (6) 可 以 求 出 UR, 将 年 均 NPP 反 向 排列 按照 
上 述 公 式 求 出 0B, 将 UF .UB 以 及 y = + 上 1.96( 置 信 
IE a =0. 05) 绘制 成 折线 图 。 通 过 分 析 统 计 序 列 
UF 和 UB 可 以 进一步 分 析 年 均 NPP 在 2000 一 2015 


(6) 


年 的 变化 趋势 ,而 且 可 以 明确 突变 的 时 间 。 如 果 UF 
和 UB 这 两 曲线 出 现 交 点 , 且 交 点 在 临界 线 y = 
+1. 96( 置 信 水 平 a =0. 05) 之 间 , 则 交点 对 应 的 年 
份 表示 年 均 NPP 开始 突变 的 年 份 。 若 UF >0, 则 表 
明 该 序列 呈 上 升 趋势 ,小 于 0 则 表明 该 序列 呈 下 降 
趋势 。 同 时 ,如 果 UF 和 UB 的 值 均 超过 临界 值 , 则 
表示 变化 趋势 显著 。 

2.2.3 地 理 探测 器 分 析 方 法 “地理 探测 器 是 探测 
和 利用 空间 分 异性 的 工具 ,最 早 应 用 于 疾病 致 病因 
ZERO” 。 本 研究 采用 因子 探测 器 和 交互 探测 器 ， 
对 塔里木 胡杨 林 国 家 级 自然 保护 区 植被 NPP 空间 
分 异 驱动 因子 进行 分 析 。 其 中 ,因子 探测 器 主要 用 于 
分 析 各 影响 因子 对 NPP 的 解释 力 的 大 小 ,公式 如 下 : 


1 m 
Pyy=1-—z >, (ny; X oh,;) (7) 


Nop i=) 


式 中 :Pyj 表 示 影 响 因 子 D 对 净 初 级 生产 力 (NPP) 
的 解释 力 ;n 为 整个 研究 区 域 样本 量 ; oz 为 整个 
区 域 NPP 的 离散 方差 ;mm 为 子 区 域内 样本 个 数 
(6=1,2,3,-+,n) 30%, ;为 子 区 域内 i 植被 NPP 的 离 
散 方差 。P ,的 取 值 范围 为 [0,1],P 值 越 大 说 明 影 
响 因 子 D 对 植被 NPP 的 解释 力 越 强 。 

交互 探测 器 主要 用 于 探测 多 个 影响 因子 交互 作 
用 后 对 植被 NPP 的 解释 力 , 即 影响 因子 X AX, 交 
互 作用 后 是 否 会 增强 或 减弱 对 植被 NPP 的 解释 力 ， 
主要 类 型 如 表 1 所 示 : 


表 1 交互 作用 探测 类 型 


Tab.1 Type of interactive detection 


判断 依据 交互 作用 
P(X, NX,) <Min(P (X,),P(X,)) 非 线 性 减弱 
Min(P (X,), P (Xs)) <P(X NX) < 单 因子 非 线 
Max(P (X,), P(X, )) 性 减弱 
P(X, NX) > Max(P (X,), P(X,)) 双 因 子 增强 
P(X NX.) = P(X,) + P(X) 独立 
P(X, AX) > P(X) + P(X) 非 线性 增强 


3 BRAM 


3.1 精度 验证 

对 区 域 或 全 球 尺 度 的 NPP 估算 结果 进行 精度 
验证 是 保证 模拟 结果 准确 与 否 的 前 提 工 作 之 一 。 通 
常 采用 两 种 方法 进行 验证 :一 是 与 实测 数据 进行 对 
比 , 二 是 将 估算 的 NPP 空间 分 布 与 他 人 在 相同 区 域 
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采用 其 他 模型 计算 的 结果 或 遥感 数据 产品 进行 对 
LE!  。 受 胡杨 林 保 护 政策 等 影响 ,不宜 对 保护 区 胡 
杨 林 进 行 破坏 ,获取 保护 区 实测 生物 量 数据 比较 困 
难 ,因此 通过 参考 单 玲 ”| . 孜 来 汗 . HEY DL 
KERSE -ER 的 硕士 学 位 论文 中 的 实测 生 
物 量 数据 进行 精度 验证 ,其 中 木 巴 热 克 … 阿 尤 普 的 
生物 量 数据 按照 胡杨 林 生物 量 模型 计算 得 到 , 取 
碳 转换 系数 0.475771 ,得 到 NPP 实测 值 。 将 NPP 实 
测 值 与 模拟 值 一 一 对 应 进行 精度 验证 (图 2) ,NPP 
实测 值 与 模拟 值 的 相关 系数 为 0.92(P < 0.01), 二 
者 线性 拟 合 曲线 斜率 为 0. 453 6, 尼 =0. 846 6(P < 
0.001) ,表明 模拟 值 与 实测 值 基本 吻合 。 该 结果 表 
明 本 研究 模拟 值 具 有 一 定 的 可 靠 性 。 

将 本 研究 NPP 佑 算 结果 与 前 人 人 研究 结果 对 比 


100 -上 
和 f y=0.453 6x+6.651 1 
R?=0.846 6 
70 


NPP 模 拟 值 / gC +m? + a 
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如 表 2 所 示 , 本 研究 2000—2015 年 年 均 植 被 NPP 
在 变化 范围 上 略 低 于 焦 伟 等 的 研究 结果 , 咯 高 于 李 
军 等 的 研究 结果 ; 焦 伟 等 发 现 2000 一 2014 年 西北 干 
旱 区 植被 NPP 呈 波 动 增长 的 趋势 ,在 2003 年 .2005 
年 .2007 年 .2013 年 出 现 大 的 峰值 ,在 2000 年 和 
2007 年 出 现 谷 值 ,与 本 研究 研究 结果 一 致 ;整体 上 
本 研究 的 模拟 值 与 不 同学 者 采用 不 同方 法 计算 得 出 
的 结果 较为 一 致 ,该 CASA 模型 可 以 很 好 的 反映 出 
塔里木 胡杨 林 国 家 级 自然 保护 区 的 植被 NPP 实际 
3.2 植被 NPP 的 空间 分 布 

2000—2015 年 塔里木 胡杨 林 国 家 级 自然 保护 
区 多 年 年 均 植 被 NPP 变化 范围 为 5. 16 ~ 303. 87 
eC .ma ,16a 平 均值 为 32.25 gC +m? .a '。 
多 年 平均 植被 NPP 以 塔里木 河 干流 为 带 向 周边 波 
动 递减 (图 3)。 


R2 同类 研究 结果 对 比 
Tab.2 Comparison of total NPP in this study 


and previous studies 


研究 时 间 模型 NPP 变化 区 间 ”文献 

2000 一 2015 CASA 5.16 ~303.87 ARM 

1981—2000 NPP-GPEM 0 ~250 [27] 

0 50 100 150 2000—2014 CASA 0 ~400 [14] 

NPP 实 测 值 / gC. m? + a? 2001—2014 MuSyQ-NPP 0 ~300 [16] 

sie ; m 1981—2008 CASA 0 ~350 28 
图 2 NPP 模拟 值 与 实测 值 的 比较 Soe 
| | i | 1982 一 2015 CASA 0 ~200 [29] 
Fig.2 Comparison of simulation results and measured NPP 
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3 2000—2015 年 年 均 NPP 空间 分 布 
Fig.3 Spatial distribution of NPP from 2000 to 2015 
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通过 分 析 塔 里 木 胡 杨 林 国家 级 自然 保护 区 不 同 a iiA a 0 
植被 类 型 (除去 无 植被 地 带 ) 年 均 植被 VPP 变化 曲 人 
线 ( 图 4) 可 知 :(1) 整体 上 看 ,保护 区 内 不 同 植被 类 90 
型 的 植被 NPP 均 呈现 增加 趋势 ,草包 与 荒漠 混交 区 p 
域 以 及 阔 叶 林 区 域 呈 现 出 明显 的 波动 增加 趋势 ,其 0 
他 植被 类 型 呈现 出 缓慢 增长 趋势 ,这 反映 出 草 旬 和 名 5 
阔 叶 林 对 干旱 区 环境 变化 更 加 敏感 。(2) 各 种 植被 六 ag 
类 型 的 植被 NPP 变化 趋势 相近 ,在 2003 年 和 2013 年 < 20| 
均 出 现 峰 值 , 在 2009 年 均 出 现 谷 值 。(3) 各 种 植被 二 
类 型 中 , 阔 叶 林 NPP 在 42. 10 ~69.56 gC +m? +a’! 
间 波 动 增长 ,年 增长 率 为 1.22 gC sm? + a! EA fe i 
NPP 变化 范围 在 25. 15 ~ 45.08 gC -m° + a 14 2000—2015 年 不 同 植被 类 型 NPP 年 际 变化 趋势 
间 ,增长 率 为 0.84 gC m” + a REAMER NPP Fig.4 Interannual changes of NPP in different 
均 较 小 ,变化 范围 均 在 20 ~35 gC.m”. a 之 间 ， biomes from 2000 to 2015 
其 增长 率 分 别 为 0.06 gC m° .a "和 0.11 
gC .marls 50 p 
3.3 MMB NPP 的 年 际 变化 | 人 

2000—2015 年 ,塔里木 胡杨 林 国 家 级 自然 保护 o 
区 植被 NPP 呈现 出 波动 增加 的 变化 趋势 (图 5) ,年 | 
均 增 加 值 为 0.523 8 gC . m+ a`, F 2003 年、 30 | 
2007 年 .2013 年 出 现 大 的 峰值 ,于 2000 年 .2009 年 25 
出 现 大 的 谷 值 ;最 高 值 出 现在 2013 年 ,其 值 为 38.9 20 Lo 
gC +m? .a ;最 低 值 出 现在 2000 年 ,其 值 为 Se ee ee 
24.55 gC +m? +a! 年 份 


利用 一 元 线性 回归 方法 得 到 塔里木 胡杨 林 国 家 
级 自然 保护 区 年 均 植 被 NPP 变 化 率 (图 6) 。 分 析 


5 2000—2015 44 


FE 均 NPP 变化 趋势 


Fig.5 Annual variation of average NPP from 2000 to 2015 
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Fig.6 NPP changing rates from 2000 to 2015 


赵 俊 红 等 :2000 一 2015 年 塔里木 胡杨 林 国 家 级 


发 现 ,2000 一 2015 年 植被 NPP 变化 率 在 -15% ~ 
15% 的 区 域 占 保 护 区 总 面积 的 40. 78% , > 15% 的 
区 域 占 保护 区 总 面积 的 16. 12% , < -15% 的 区 域 
占 保护 区 总 面积 的 43.1% ,表明 16 a 间 保 护 区 植被 
NPP 变化 较为 稳定 。 塔 里 木 河 干流 的 植被 NPP 呈 
明显 上 升 趋势 ,大 部 分 地 区 NPP 变化 率 >30% , E 
要 植被 有 阔 叶 林 和 草 匈 ;面积 低 于 -30% 的 区 域 主 
要 集中 在 保护 区 的 东北 部 ,主要 植被 类 型 为 荒漠 、 
草 甸 。 

将 塔里木 胡杨 林 国 家 级 自然 保护 区 2000 一 
2015 年 年 均 植 被 NPP 值 进 行 Mann-Kendall 趋势 分 
析 和 突变 检验 ,算出 统计 量 Z =2. 29 >0( 置 信 水 平 
a =0. 05) ,表明 塔里木 胡杨 林 国 家 级 自然 保护 区 年 
均 植 被 NPP 呈现 增加 趋势 ,图 7 中 UF 统计 量 代 表 
植被 NPP 的 Mann-Kendall 顺序 统计 曲线 ,UB 统计 
量 代 表 植 被 NPP 的 Mann-Kendall 道 序 统计 曲线 ， 
0.05 显 著 水 平 为 临界 线 。 由 UF 和 UB 曲线 可 知 : 塔 
里 木 胡 杨 林 国家 级 自然 保护 区 2000—2015 年 年 均 
植被 NPP 的 UF >0 ,表明 塔里木 胡杨 林 国 家 级 自然 
保护 区 年 均 植 被 NPP 呈 增 加 的 趋势 ;UF 和 UB hh 
线 于 2001 一 2002 年 和 2011 一 2012 年 出 现在 临界 线 
之 间 有 两 个 明显 的 交点 ,表明 胡杨 林 自 然 保 护 区 年 
均 植被 NPP 在 2001 一 2002 年 和 2011 一 2012 年 均 出 
现 突 变性 上 升 ,分 别 增加 了 4.46 gC -m° + a! Al 
3.15 gC +m +a_';UF 和 UB 的 值 于 2007—2008 
年 均 超 过 临界 值 ,表明 塔里木 胡杨 林 年 均 植 被 NPP 
在 2007 一 2008 年 有 突变 特征 ,突变 年 份 之 间 的 年 均 
植被 NPP 下 降 了 5.96 gC +m? .a 。 
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图 7 2000—2015 年 年 均 NPP 变化 率 
Fig.7 NPP changing rates from 2000 to 2015 
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3.4 植被 NPP 空间 差异 影响 因素 分 析 

气候 变化 特别 是 降水 和 气温 的 变化 ,直接 影响 
植物 光合 作用 及 新 陈 代谢 而 改变 植被 生产 力 [9 。 
不 同 土地 利用 /覆盖 类 型 对 应 不 同 的 土地 生态 系统 
结构 ,具有 不 同 的 植物 生产 功能 "1 。 不 同 土壤 类 型 
在 颗粒 状 矿 物质 ` 有 机 物质 ,水 分 .空气 和 微生物 等 
的 不 同 ,所 产生 的 土壤 肥力 不 同 ,直接 影响 植被 生长 
PREC 。 同 时 ,在 塔里木 河 干流 沿岸 ,径流 是 限制 
胡杨 生长 的 主要 因素 "””。 据 此 ,选取 2000—2015 
年 多 年 平均 NPP 作为 被 探测 因子 ,选取 气候 因子 
(包括 降水 .气温 和 莹 散发 因子 ) .土地 利用 /覆盖 、 
河流 缓冲 区 因子 (距离 河道 的 距离 ) .土壤 类 型 作为 
解释 因子 。 

对 气候 因子 进行 自然 间断 点 分 级 聚 类 ,对 土地 
利用 /覆盖 以 及 土壤 类 型 的 分 类 采用 相应 的 分 类 标 
准 ,河流 因子 根据 缓冲 区 进行 分 类 ,采用 地 理 探测 器 
分 析 方 法 ,得 到 各 个 解释 因子 对 多 年 平均 NPP 动态 
变化 的 决定 力 尸 值 及 其 显著 性 水 平 。 

由 表 3 可 知 :所 有 变量 都 通过 了 0.01 水 平 显 著 
性 检验 ,表明 上 述 因 素 对 NPP 均 有 显著 的 决定 力 。 
从 决定 力 因子 的 大 小 来 看 :(1) 土地 利用 /覆盖 对 植 
被 NPP 变化 影响 作用 最 大 ,决定 力 达 0.5。 土 地 利 
用 /覆盖 类 型 及 面积 决定 了 生态 系统 结构 和 功能 ,后 
时 ,土地 利用 /覆盖 是 决定 该 区 域 蒸 散发 的 重要 因 
K ,对 植被 NPP 有 重要 影响 。(2) 在 气候 因子 
中 , 莹 散发 及 降水 对 植被 NPP 变化 的 决定 力 分 别 为 
0.22 和 0. 19 ;而 气温 对 植被 NPP 变化 的 决定 性 较 
小 ,只 有 0.07。 蒸 散发 是 影响 区 域 水 循环 过 程 的 重 
要 因子 之 一 ,近年 来 该 区 域 蒙 散发 出 现 增 大 趋 
ADU ,直接 影响 了 区 域 水 资源 时 空 配置 与 平衡 。 保 
护 区 近 82.61% 的 地 区 与 降水 呈现 正 相 关 关 系 , 植 
被 NPP 与 降水 呈 高 度 正 相关 (0.5 < R<1) 和 中 度 
正 相关 (0.3 <R<0.5) 的 区 域 的 面积 约 占 整个 保护 
区 面积 的 49.21% 。 保 护 区 约 63. 02% 的 区 域 与 气 
温 呈 现 弱 相关 关系 ,这 可 能 是 由 于 温度 变化 对 胡杨 
生长 的 影响 具有 “滞后 效应 ”, 年 平均 温度 的 变化 可 


表 3 各 影响 因素 对 年 均 NPP 的 决定 力 地 理 探测 结果 
Tab.3 Interactive impact of factors on average NPP 


with geographic detector analysis 


ek an eww 十 地 利用 河流 缓 CE 


决定 力 0.19 0.07 0.22 0.50 0.18 0.09 
显著 性 水 平 0.00 0.00 0.00 0.00 0.00 0.00 


tó wy 


AES F — AF oe aE 。(3 ) 
河流 缓冲 区 因子 对 植被 NPP 变化 影响 较 高 ,决定 力 


KA, KNA HAI (0.56) RBE N EIA 
(0.56) 以 及 土地 利用 站 土壤 类 型 (0.56) 交 互 后 


为 0.18。 塔 里 木 河 干流 荒漠 河岸 林 的 重要 保护 范 
围 在 0.8 ~12.3 km’? ,结合 NPP 空间 分 布 , 对 河流 
做 8 km 缓冲 区 ,对 时 令 河 做 2 km 缓冲 区 ,把 河流 组 
冲 区 与 多 年 平均 植被 NPP 进行 空间 统计 , 发现 河流 
缓冲 区 内 植被 NPP 总 量 占 保护 区 植被 NPP 总 量 的 
88.08% ,平均 值 为 45.06 gC m° a !, 约 为 全 区 
植被 NPP 均值 的 1.4 倍 ;其 中 ,塔里木 河 干流 植被 
NPP 总 量 占 保护 区 植被 NPP 总 量 的 74. 16% ,平均 
值 为 56. 20 gC m? a! , 约 为 全 区 年 均 NPP 的 
1.7 o (4) 土壤 类 型 对 植被 NPP 变化 影响 较 小 ， 
这 是 由 于 保护 区 土壤 颗粒 以 粉 粒 和 砂粒 为 主 , 胡 杨 
在 粉 土 及 沙土 中 发 育 不 良 甚至 衰败 ;土壤 有 机 质 含 
量 最 高 为 0.8% ,土壤 养分 贫 将 ;同时 ,土壤 盐 碱 化 
是 抑制 胡杨 萌发 与 生长 的 主要 因素 扩 ; 。 
植被 NPP 变化 是 多 种 因素 相互 作用 的 结 

因此 需要 进一步 研究 多 种 因素 交互 作用 对 植被 
NPP 变化 的 影响 。 利 用 地 理 探测 器 进行 了 交互 作 
用 探测 ,探测 结果 ( 表 4) 表 明 : 从 交互 作用 的 效果 
来 看 ,保护 区 植被 NPP 变化 影响 因素 交互 作用 主 
要 表现 为 非 线性 增强 和 双 因 子 增强 ,这 说 明 任意 
两 个 因素 交互 后 对 植被 NPP 变化 的 决定 力 均 会 显 
著 提 升 ,也 就 是 说 保护 区 植被 VPP 的 变化 受到 降 
水 .气温 . 蒸 散 发 .土地 利用 /覆盖 .河流 缓冲 区 、 土 
壤 类 型 等 因子 的 共同 制约 ,两 个 因素 交互 作用 后 
的 决定 力 明 显 高 于 原来 单 因 素 ; 从 交互 作用 大 小 


因子 决定 力 最 大 ,其 次 为 气温 站 燕 散发 (0.55 ) ,这 
表明 它们 的 交互 作用 相对 于 其 他 因子 之 间 的 交互 
作用 较为 显著 。 

在 干旱 区 ,地 下 水 对 于 植物 生长 起 着 至 关 重 要 
的 作用 ,地 下 水 位 埋 深 影响 着 植被 NPP 的 大 小 。 
此 ,本 研究 利用 徐 梦 尾 等 ”在 塔 河中 游 监 测 断 面 实 
测 数 据 , 分 析 2001 一 2012 年 地 下 水 位 埋 深 与 植被 
NPP 的 关系 。 从 图 8 看 出 , 除 2007 年 以 及 2009 一 
2012 年 外 ,其 余年 份 地 下 水 埋 深 与 年 均 植被 NPP 
总 体 上 呈 相 反 的 变化 趋势 , 即 随 着 地 下 水 位 埋 深 的 
增加 ,年 均 植被 NPP 不 断 下 降 ,2003 年 地 下 水 位 埋 
深 最 浅 , 年 均 NPP 达到 最 大 值 ;2001 一 2003 年 ,地 
下 水 埋 深 保持 在 3 m 左右 ,年 均 植 被 NPP 在 26 ~ 
35 gC +m + a ' 间 波动 ;2003 一 2009 年 ,地 下 水 埋 
深 持续 增加 ,年 均 植 被 NPP 从 35 gC sm? +a! F 
降 到 26 gC sm? .a '。 总 体 来 说 ,地 下 水 埋 深 的 
时 间 变 化 影响 着 植被 NPP 的 年 际 动 态 变 化 。 
3.5 讨论 

胡杨 林 作 为 新 疆 干 旱 荒 漠 区 河岸 林 的 关键 树 
种 ,对 维持 荒漠 生态 系统 具有 十 分 重要 的 作用 。 近 
年 来 , 受 生态 输 水 的 影响 , 塔 河 下 游 胡 杨 林 不 断 恢 
复 。 相 关 研 究 表明 ,2016 年 塔里木 河 下 游 的 英 苏 、 
喀 尔 达 依 、 阿 拉 干 以 及 依 干 不 及 麻 断 面 生 物 量 分 别 
为 626.306 kg . hm ,208. 632 kg . hm ,555. 085 
kg + hm- ,276.428 kg . hm?) ,折算 成 NPP 分 


表 4 植被 NPP 影响 因子 交互 作用 
Tab.4 Interactive impact of factors on NPP 


X NX, P(X NX,) 判断 依据 交互 作 
降水 (0. 19) 气温 (0.07) 0. 46 P(X, NX,) >P(X,) +P(X3) 非 线性 增强 
降水 (0. 19) ARAL (0. 22) 0. 43 P(X NX,) >P(X,) + P(X, ) 非 线性 增强 
降水 (0. 19) A 土地 利用 /覆盖 (0. 50) 0. 56 P(X, OX) > Max( P(X, ) ,P(X,)) 双 因 子 增强 
降水 (0. 19) 站 河流 缓冲 区 (0. 18) 0.29 P(XINX,) > Max(P(X,) ,P(X,)) 双 因 子 增强 
降水 (0.19) 路 土壤 类 型 (0. 09) 0. 35 P(X, NX,) >P(X,) + P(X) 非 线性 增强 
气温 (0.07) NZR (0. 22) 0. 55 P(X, OX) >P(X,) + P(X) 非 线 性 增强 
气温 (0.07) A 土地 利用 /覆盖 (0. 50) 0. 57 PEna)=P(X)+P(X2 ) 独立 
气温 (0.07) 河流 缓冲 区 (0. 18 ) 0.31 P(X, NX,) >P(X,) +P(X3) 非 线性 增强 
气温 (0.07) N 土壤 类 型 (0. 09) 0. 22 P(X, NX,) >P(X,) + P(X, ) 非 线性 增强 
ARB (0. 22) mm 土地 利用 /覆盖 (0.50) 0. 56 P(XINX,) >Max(P(X,) ,P(X,)) 双 因 子 增强 
FRB (0. 22) Nn 河流 缓冲 区 (0. 18 ) 0. 29 P(X, 1X, ) >Max(P(X,) ,P(X,)) 双 因 子 增强 
REE (0. 22) N 土壤 类 型 (0. 09) 0. 40 P(X, 0X5) >P(X,) +P(X,) 非 线性 增强 
土地 利用 /覆盖 (0. 50) 站 河流 缓冲 区 (0. 18 ) 0. 52 P(X 1X, ) >Max(P(X,) ,P(X2)) 双 因 子 增强 
土地 利用 /覆盖 (0. 50) m 土壤 类 型 (0. 09 ) 0. 56 P(X, 1X, ) >Max(P(X,) ,PCD )) 双 因 子 增强 
河流 缓冲 区 (0. 18) N 土壤 类 型 (0.09 ) 0. 26 P(X, NX, ) > Max( P(X, ) ,PCX2 ) ) 双 因 子 增强 
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图 8 2001 一 2012 年 地 下 水 埋 深 与 植被 NPP 的 关系 
Fig.8 Relationship between groundwater depth and vegetation 
from 2001—2012 


别 为 29.75 gC - m`” »a`™.9.91 gC: m” .a !、 
26.37 gC +m? +a! 13.13 gC e-m” .a ,塔里木 
河 下 游 NPP 均值 为 19.79 gC m° a~ RPK 
2000—2015 年 多 年 年 均 NPP H 32.25 gC +m” ， 
a! , 均 高 于 下 游 4 个 断面 的 NPP, 这 反映 生态 输 水 
起 到 了 积极 作用 ,保护 区 的 保护 效果 整体 上 较为 良 
好 。 

生态 系统 中 不 同 组 分 植被 NPP 可 以 反映 其 在 
生态 系统 中 的 重要 性 ,我 国 不 同 森 林 类 型 植被 NPP 
也 存在 差异 ,我国 落叶 针 叶 林 、 常 绿 针 叶 林 、 落 叶 阅 
叶 林 年 均 NPP 分 别 为 438.8 gC +m? + a7! 367.1 
eC: m+ a! 642.9 eC: m+ a!U] ,新 疆 天 山 
以 及 阿尔 泰山 针 叶 林 NPP 在 200 ~400 gC .m 一 + 
a 105。 研究 区 年 均 NPP 为 32.25 gC +m? +a”, 
变化 范围 在 5. 16 ~ 303.87 gC +m? +a 之 间 , 均 
低 于 落叶 针 叶 林 、 常 绿 针 叶 林 、 落 叶 阔 叶 林 等 植被 类 
型 ,也 低 于 同属 于 干旱 区 的 天 山 和 阿尔 泰山 ,这 说 明 
气候 条 件 特 别 是 降水 是 塔里木 河 植 被 生长 发 育 的 主 
要 限制 条 件 。 黑 河流 域 2008 一 2014 年 植被 NPP 变 
化 范围 在 0 ~ 560.5 gC + m°? + a~! ZA, EES 
于 42.25 ~46.83 gC > m? + a!) ,研究 区 植被 
NPP 演化 状况 与 黑河 流域 较为 接近 。 因 此 ,在 干旱 
区 脆弱 生态 系统 中 ,胡杨林 对 于 森林 生态 系统 物质 
与 能 量 循环 过 程 以 及 碳 汇 具 有 重要 意义 。 


4 结论 


从 保护 区 尺度 出 发 ,基于 植被 VPP 气温 、 降 
水 .土地 利用 /覆盖 河流 等 数据 ,采用 趋势 分 析 、 
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P3 然 保 护 区 NPP 时 空 动态 特征 及 其 影响 因素 


Mann-Kendall .相关 分 析 、 地 理 探 测 器 等 方法 ,分 析 
了 塔里木 胡杨 林 国 家 级 自然 保护 区 植被 NPP 时 空 
分 布 及 其 差异 ,探测 了 影响 植被 NPP 空间 分 异 的 因 
素 ,结果 表明 : 

(1) 2000 一 2015 年 塔里木 胡杨 林 国 家 级 自然 
保护 区 多 年 年 均 植被 NPP 在 5.16 ~ 303. 87 
gCom”? +a! 之 间 呈 波动 增长 趋势 ,平均 值 32. 25 
“a”, 年 均值 在 2013 年 出 现 最 大 值 ,在 
2000 年 出 现 最 小 值 。 空 间 上 ,年 均 NPP 呈现 出 以 
塔里木 河 干流 为 带 向 周边 波动 递减 的 空间 分 布 特 
征 。 不 同 植被 类 型 的 年 均 NPP 均 呈 增加 趋势 ,在 
2003 年 和 2013 年 出 现 峰值 ,在 2009 年 出 现 谷 值 ; 
各 种 植被 类 型 中 , 阔 叶 林 植被 NPP 最 高 , 草 多 次 
之 , 灌 从 和 荒漠 NPP 均 较 小 ;就 增长 速度 而 言 , 阔 
叶 林 年 均 NPP 增长 速度 最 快 ,其 次 为 草 旬 , 灌 从 最 
慢 。 

(2) 保护 区 近 16 a 植被 NPP 呈现 波动 增长 趋 
势 ,年 均 增加 值 为 0.523 8 gC .mm 一. a ', 2000— 
2015 年 间 塔 里 木 胡杨 林 植 被 NPP 出 现 突变 特征 ， 
在 2001 一 2002 年 和 2011—2012 年 均 出 现 突变 性 上 
升 , 在 2007 一 2008 年 有 突变 性 下 降 。 

(3) 决定 保护 区 植被 NPP 空间 差异 的 各 因子 
重要 程度 依次 为 土地 利用 /覆盖 、 蒸 散发 、 降 水 、 河 
流 .土壤 类 型 .气温 ,说 明 土 地 利用 / 履 被 变化 对 植被 
NPP 变化 影响 最 大 ,其 次 水 资源 条 件 是 干旱 区 生态 
系统 的 关键 因子 ,气温 条 件 影响 最 小 ;同时 ,地 下 水 
埋 深 影响 着 植被 VPP 的 年 际 动态 变化 。 任 何 两 个 
因子 交互 作用 后 ,对 植被 NPP 的 解释 力 均 大 于 单一 
因子 的 解释 力 ,说 明 多 因子 协同 作用 造成 塔里木 胡 
杨 林 植 被 NPP 的 空间 分 异 。 


eC +m 


5 建议 


根据 近 16 a 保护 区 植被 NPP 演化 状况 及 其 制 
约 因素 ,对 其 生态 恢复 与 保护 提出 以 下 几 个 建议 : 

(1) 加 强 塔 河流 域 综合 治理 ,优化 塔 河上 中、 
下 游 水 资源 配置 ,协调 下 游 生态 输 水 和 中 游 生态 需 
水 。 

(2) 加 强 保 护 区 河流 两 侧 地 下 水 位 动态 监测 ， 
RRES VA YO PEK BF WE, FPA 
河道 古河 汉 ,实施 引 潜 灌溉 ,促进 保护 区 胡杨 林 等 
自然 植被 的 恢复 和 更 新 。 
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(3) 对 于 距离 河道 较 远 以 及 防护 堤 外 的 区 域 ， 


可 种 植 人 工 幼 苗 , 以 弥补 自然 成 苗 困 难 的 局 面 。 
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Temporal-spatial characteristics and influencing factors of the vegetation 
net primary production in the National Nature Reserve of Populus 
euphratica in Tarim from 2000 to 2015 


ZHAO Jun-hong'*, ZHOU Hua-rong', LU Ya-yan'?, SUN Qing-xiang'” 
(1 Xinjiang Institute of Ecology and Geography , CAS , Urumqi 830011 , Xinjiang, China; 
2 University of Chinese Academy of Science , Beijing 100049 , China ) 


Abstract; Net Primary Productivity (NPP) is an important indicator to characterize the production capacity and 
quality of terrestrial ecosystems. Previous studies are mainly based on large and middle scale units and pay more at- 
tention to semi-humid and humid areas, therefore we select National Nature Reserve of Populus euphratica in Tarim 
Basin , Xinjiang ,China as the study area. Based on the NPP , climate , land use/land cover ,river data and using trend 
analysis , Mann-Kendall test, correlation analysis and geographical detector methods, this paper analyzed the tempo- 
ral-spatial characteristics and influencing factors of the vegetation NPP in the National Nature Reserve of Populus 
euphratica in Tarim from 2000 to 2015. The results are as follows: First, the average vegetation NPP from 2000 to 
2015 was 32.25 gC - m“? + a’! and ranged from 5. 16 —303. 87 gC - m~? + a~’. The average NPP was declined as 
the belt form from the two wings of Tarim River to the place far away from the river. Different vegetation types con- 
tribute differently to annual vegetation NPP growth, and NPP tended to increase the most for broadleaf forest , but 
least for shrublands. Second ,the NPP followed an overall trend of fluctuating growth with 0.523 8 gC +m? + a7! 
every year. The NPP had a sharp growth trend during the period from 2001 to 2002 and the period from 2011 to 
2012 and drastically reduced during the period from 2007 to 2008. Last, there are many factors affecting NPP , main- 
ly including land use/land cover , evapotranspiration , precipitation and river buffer zone ,and the NPP changes were 
mainly caused by the synergistic effect of multiple factors. The groundwater depth had a significant influence on veg- 
etation NPP. This study provides decision support for the ecological environment protection and governance of Na- 
tional Nature Reserve of Populus euphratica in Tarim Basin, and provides a reference frame for maintaining the eco- 
system sustainability of the green ecological corridor in Tarim River. 
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Populus euphratica in Tarim 


